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Generation of wave structures by a two-dimensional object (laser beam) moving in a two- 
dimensional two-component Bose-Einstein condensate with a velocity greater than both sound ve- 
locities of the mixture is studied by means of analytical methods and systematic simulations of the 
coupled Gross-Pitaevskii equations. The wave pattern features three regions separated by two Mach 
cones. Two branches of linear patterns similar to the so-called "ship waves" are located outside the 
corresponding Mach cones, and oblique dark solitons are found inside the wider cone. An analytical 
theory is developed for the linear patterns. A particular dark-soliton solution is also obtained, its 
stability is investigated, and two unstable modes of transverse perturbations are identified. It is 
shown that, for a sufficiently large flow velocity, this instability has a convective character in the 
reference frame attached to the moving body, which makes the dark soliton effectively stable. The 
analytical findings are corroborated by numerical simulations. 

PACS numbers: 03.75.Kk 



I. INTRODUCTION 



Breakdown of superfluidity at large velocities of the flow is a fundamentafly important problem which has been 
widely studied in physics of quantum liquids, such as ^He and Bose-Einstein condensates (BECs) of dilute atomic 
gases (see book Jj and references therein). As is known, the breakdown is caused by opening of channels of radiation 
of elementary excitations in the fluid, as it happens, for instance, with increase of the velocity of a body ("obstacle") 
moving through the superfluid (in the experiment, a far-blue-detuned laser beam, which produces a local repulsive 
force acting on atoms, usually plays the role of the obstacle 0,3). In the BEG, solitary waves and vortices are readily 
generated if the size of the obstacle is of the order of (or greater than) the characteristic healing length of condensate. 
The generation of these structures manifests itself as an effective dissipation, which implies the loss of superfluidity at 
some critical value of the obstacle's velocity On the other hand, if the size of the obstacle is much smaller than 

the healing length, the main loss channel corresponds to the Cherenkov radiation of Bogoliubov excitations, and it 
opens at supersonic velocities of the obstacle Q . If a large obstacle moves at a supersonic velocity, then the amplitude 
of the generated waves becomes large too. In the latter case, two dispersive shocks, which start their propagation 
from the front and the rear parts of the moving body, are formed. Far from the body, the shock front gradually 
transforms into a linear "ship wave" located outside the Mach cone 0, H, lipt - whereas the rear zone of the shock is 
converted into a "fan" of oblique dark solitons located inside the Mach cone llll,[l3,[l|. Although, as is well known, 
such dark solitons are unstable with respect to transverse perturbations, it was shown in that for a flow velocity 
greater than some critical value this instability becomes a convective one (rather than being absolute) in the reference 
frame moving along with the obstacle. This fact actually means that the dark solitons are effectively stable in the 
region around the obstacle. Some of these structures have already been observed in experiments [ A Il5 | . and similar 
non-stationary dispersive shocks were studied both theoretically and experimentally in BEG [3. [l6l.]l7l fisl [T9| (see 
also review [Ifll) and nonlinear optics [HI, [H, [H [H, HI [Hi . 

The picture described above corresponds to a single-component BEG whose mean-field dynamics is governed by the 
respective Gross-Pitaevskii equation [l(]. At the same time, one of the focal points of activity in B EG phys ics has been 
the study of multi-component settings, a prototypical one b eing presented by binary mixtures |32l. [ssl. [33 |. These 
two-component media exhibit phase-separation phenomena [35l [sgI . \3T\ due to the nonlinear interaction between 
the different atomic species (or possibly different hyperfine states) that constitute the mixture. The formation of 
robust single- and multi-ring patterns [ssl Issj , the evolution of initially coincident triangular vortex lattices through 
a turbulent regime into an interlaced square vortex lattice [s^ in coupled hyperfine states in the ®^Rb condensate, 
and the study of the interplay between atomic states at different Zeeman levels in the ^"^Na condensate, forming 
striated magnetic domains in optical traps [dol . l4l| , are only a small subset among the many possibilities that multi- 
component BEGs can offer. It should also be noted that mixtures with a higher number of components, namely spinor 
condensates [12, are known too. They have been realized with the help of far-off-resonant optical techniques for 
trapping ultracold atomic gases |4^] , which, in turn, allowed the spin degree of freedom to be explored (previously, it 
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was frozen in magnetic traps). 

Our aim in the present work is to unite these two areas by investigating the effects generated by the motion of an 
obstacle in a "pancake" -shaped, i.e., effectively two-dimensional (2D), two-component BEG. In an earlier work [is! ]. 
the critical situation, when the obstacle had the velocity comparable to the two speeds of sound in the two-component 
system, was examined. Here, we extend the analysis to the supercritical case, when the speed of the moving body may 
be significantly higher than the sound speed(s). We demonstrate that two branches of linear, so-called "ship- wave", 
patterns form outside of the two Mach cones (which are associated with the speeds of sound), while oblique dark 
solitons are located inside the wider Mach cone. While these dark solitons are unstable at relatively low velocities of 
the motion of the obstacle, they can be convectively stabilized at sufficiently high values of the velocity. 

The presentation is structured as follows. In section II, we outline the model under consideration. In section III, 
we study the linear (ship) waves by means of both analytically and numerical methods. In section IV, we study dark 
solitons (again, by dint of both analytical and numerical methods). Finally, we summarize our findings in section V, 
and discuss potential directions for future work. 



II. MATHEMATICAL MODEL AND SETUP 



In the usual mean-field approximation, the 2D flow of a binary condensate mixture past an obstacle obeys the 
system of nonlinearly coupled Gross-Pitaevskii eq uations (GPEs), with spatial coordinates r — {x,y). In the scaled 
form, the equations take a well-known form P, [2^ , 

= -^At/^i + (511 + 512 |V^2|')^i + Vir, i)V'i, 
dt 2 

dip2 1 2 2 

i-Q^ = -^^1p2 + {9l2\lp2\ +922\lp2\ )'02 + V(r, i)V'2- 

We assume that atoms in both species have the same mass (i.e., they represent different hyperfine states of the same 
atom, as is often the case, see e.g., Ref. [S^ and references therein), V{r, t) being the potential induced by the moving 
obstacle, which is identical for both components. 

Small-amplitude waves and solitons correspond to potential fiows with zero vorticity, for which Eqs. ([T]) can be 
transformed into a hydrodynamic form by means of substitutions 



ipi (r, t)^ (r, t) exp ^ Ui {r' ,t)dr' - ifiit^ , ^2(1",^) = \/n2(r,t)exp J U2(r', t)dr' - i/i2i^ , 



(2) 



where ni^2(r, i) are the atom densities of the two BEG components, Ui_2(r,t) their velocity fields, and /ii_2 the 
respective chemical potentials. Substituting Eqs. ([2]) into the Eqs. IJ) we arrive at the following system, 

ui^t + V • (niui) 0, 7^2, t + V • (n2U2) = 0, 
ui,t + (ui • V)ui + giiVni + gi2^n2 + V f - ^\ + Vl/(r, t) = 0, 



8nf 4ni y ^ ' ' ' (3) 

U2,t + (U2 • V)u2 + gi2Vni -I- 522 Vn2 -I- V {^^^ " + '^^('^' ^) " 0. 

The first pair of the equations represents the conservation of the number of atoms in each component, and the second 
pair corresponds to the Euler's equations for fiuid velocities under the action of the pressure induced by interactions 
between atoms, the obstacle's potential, and quantum dispersion. 

We consider waves generated by the obstacle moving at constant velocity U along the x-axis, 

V{v,t)^V[x-Ut), (4) 

through a uniform condensate, so that at \x\ 00 both components have constant densities and vanishing velocities. 
This setting implies the absence of a trapping potential in the plane of the quasi-2D BEG (or, more realistically, a 
very weak trapping). In the reference frame moving along with the obstacle, the unperturbed condensate fiows at 
constant velocity u = (— [/, 0), hence the respective boundary conditions for the densities and velocities are 

ni Tiio, 712 «20, ui ^ u, U2 — > u at |x| — > 00. (5) 

As follows from Eqs. (O, chemical potentials /ii^2 are related with the asymptotic densities 7iio,?i20 by the relations 

/ll = 5ll"-10 + 3l2»l20, = 312^110 + 522?l20- (6) 



In this reference frame, the wave pattern is a stationary one, which is convenient for analytical considerations. 
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III. LINEAR "SHIP WAVES" 

We first aim to consider linear waves generated by the moving obstacle, assume that potential V{r, t) is weak 
enough to apply the perturbation theory Q based on linearized equations. Actually, the approximation is valid 
everywhere where the amplitude of the waves is small — in particular, far enough from the obstacle outside the Mach 
cones associated with the speeds of sound. Thus, we introduce small deviations from the uniform state, namely, 

ni^nio + n[, n2=n-20 + n2, Ui^u + u'^, U2 = u + U2, (7) 

and accordingly linearize Eqs. ([3]): 

n'it+nio(V-u;) + (u- V)n; = 0, ^ + n2o(V • u^) + (u • V)n^ = 0, (8) 



(9) 



u'l , + (u • V)u; + giiVn'i + g^^n'^ ~ -^V{An[) = -VV, 

U2,t + (u • V)u'2 + gi2Vn; + g22Vn^ - -^V{An'^) = -VV. 

In the stationary case, all time derivatives vanish. Further, we apply the Fourier transform, 

n[ir,t)^ J J n;(k,i)e*"-^. (10) 

Then, we eliminate the velocity perturbations in Eqs. ^ in favor of the densities, arriving at the following system of 
two linear equations: 

[-(k • u)^ + fc2(5iinio + k^/4:)]fi[ + k^gunion^ = -fc^fmo, ^^^^ 
k^gi2n2on[ + [-(k • u)^ + fc2(522"20 + fcV4)]"2 = ~k^Vn2o, 

where tildes denote the Fourier components. This linear system can be readily solved for h'l 2 and the inverse Fourier 
transform yields 



-nw 



-n.20 



k^V{[-{k ■ U)2 + fc2(ff22n20 + fcV4)] - fc2ffi2»20} d^k 

((k-u)2-cj2)((l^.y)2_^2) (27r)2^ 

fc2F {[-(k ■ u)2 + fc2(gnnio + fc2/4)] - fe2gi2»io} d^k 

((k-u)2-u;2)((k.u)2_^2) (2,^)2^ 

where the dispersion relations for the linear waves in the binary mixture are given by 



(12) 



/ >2 _ 11,2 



giiniQ + g22n2Q + \k'^ ± \l {giiuw - .922^20)^ + ^gi2i^wn2o 



(13) 



In the long-wave limit, Eqs. (fT3)) yield the sound velocities [45j], c± — uj{k)/k, given by: 



uj{k) / giinio + 522»^20 ± V {gimw - 522"-2o)^ + 4gj2"io"20 

^±"lToT = V 2 ■ 

Expressions (fT^ can be analyzed following the lines of Ref. [§, [l3|. To this end, we introduce new coordinates (see 
Fig. 1), 

r = (— r cosx, J'sinx), k = (fccosry, fcsin?]), (15) 

and assume that the wavelength of the pattern is much greater than the characteristic size of the obstacle, hence its 
potential can be approximated by the form V{r) = Vbi5(r). Then, we obtain 

, _ Wpnw r r fc[(.922 ~ 9i2)n2o + k^A - cos^ Tyje'*^-- 
.2 J ik^-kl-^OW-k^_-^0) 

, _ mn20 r r fc[(gii - gi2)nio + fcV4 - £/2 cos2 7^^^'^ 

^2 J (fc2_^^_,0)(fc2_fc2 _,0) 



FIG. 1: (Color online.) Coordinates defining radius-vector r and wave vector k. The latter one is normal to the wave front of 
one of the ship- wave modes, which is shown schematically by a curve. 



where infinitesimal imaginary parts in the denominators are written to define correct contributions from the poles 
corresponding to adiabatically slow switching on of the potential, and 



fc± = 2y'C/2cos2r;-4. (17) 

Next, we split the integration domain of ry into two parts: /^^yj ^ /-^/2 ^ I^J'z'^ '^'^ ^^er replacement 
7]' = ?7 — TT in the second term, one can notice that the integrand turns into its own complex conjugate, allowing one 
to write the integrals in Eq. (fT6| as 

, _ SFomo r fc[(.g22 ~ .912)^20 + feV4 - cos^ ^]e""- 

, _ SVpn^o ^ r'^ , r fc[(gii - .gi2)nio + ^^4 - L/^ cos^ ^Y'^^ 

-2 ^'y_//''io (fc2_fc2 „,0)(P-P _,0) 

The integration over fc should be carried out over the positive real half-axis. However, we can add to this path a 
quarter of an infinite circle and the imaginary half-axis in the complex fc plane, to build a closed integration contour. It 
is easy to show that the contribution from the quarter of the infinite circle is zero, and contribution from the imaginary 
axis depends on r as r"^, decaying at large r much faster than the contribution of the pole, which is ~ r^^l"^ . Thus, 
far from the obstacle, it is sufficient to keep only the contribution from the poles, which yields (i^ = tt — x — ?/) 



2Vonio 1 r 2 / NIT r'"^ 

" - (522 - 5i2)"2oJIm 



7r(4-c2_) 

2Vbn20 1 r 2 / ^ IT C ^'^ 



f d,7e"=+'^^°^'' - [cl - (322 - 5i2)"2o]Im / dr,e 

J-tt/2 J-tt/2 




(19) 



[4 - (511 - 5i2)"io]Im r d?7e''+'^^°^'' - [c^- - (.gii - <?i2)ni 

J-7r/2 

Far from the obstacle, where phases k±r — rs± are large, we are dealing with large values of 

s±iv) = k±{r])cos{x + v), (20) 

and the integrals in Eq. (|19p can be estimated by means of the standard stationary-phase method. Since calculations 
of both integrals are identical, we consider, for definiteness, the integral over Condition ds+fdrj = is an equation 
for the stationary-phase point, which can be easily transformed to 

tanz^+ = (2J7^/fc^) sin2r/+, (21) 

or, with regard to the definition of v, we obtain an expression for 

(l + fc^/(24))tan7y+ 
*""^+=C/V4- (1 + ^/(24))- 
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|U^ 1^ when v=3.5, t = 35 when v=3.5, t = 35 




-50 50 100 -50 50 100 

X X 



FIG. 2: Spatial contour plots of densities of the two components (left and right panels, respectively) of the BEG for the flow 
velocity U = 3.5, at t = 35. The form of the obstacle'e potential used in this case is V{x, y,t) — 2 j^sech ^\/ (x — Uty + y2/2^ j . 
Solid and dashed thin lines correspond to wave crests of two modes of ship waves. 



Here 77+ takes values in interval 



arccos < ri+ < arccos , M+ = — , (23) 

M+ - - M+ c+ 



while the corresponding vector {x,y}, given by the parametric expressions, 

x{t]+) = -^{Ml cos2tj+~1)cos7j+, y{r^+) = -^{2MIcos^ ^+ - l)sm7j+, (24) 

fV I ru j 

moves along a curve with constant phase (e.g., a crest line). 

As usual in the method of stationary phase, we reduce the integrals in (jl9p to Gaussian ones and, as a final result, 
obtain 



2Vbnio I r 2 / ^ 1 / 

44-c^_) \[^+'("^^~g"^"^"]Vfc77 

[l + (4£/Vfc^)sin^ (27?+)]'^" 

[1 + (AU^/kl) cos 27?+ + (l2C/4/fc|) sin^ (27]+)]' 

^ (l+(4[/Vfc^)sin^2ry„)'^' 



-7- cos(fc+r COSZ/+ — 7r/4) 
1/2 



(25) 



- [c_ - (.922 - 5l2)?^2oJ 



X cos(fc_r cos U- — 7r/4) 



^-^ [1 + (AU^/kl) cos (277_) + {l2U^/ki) sin^ (277.)] 



where i>± = n — x± — '?± , and similar expressions can be derived for the density oscillations of the second component 
of the binary condensate. 

An example of wave patterns generated by numerical simulations of Eqs. ([T]) is displayed in Fig. 2, for parameters 
gii — 1.5, 322 = 1.25, gi2 = 1.0, and nio = 1.0, 7120 = 2.0 and velocity U — 3.5. Three different types of waves can 
clearly be distinguished in this picture. First, one observes ship waves outside the Mach cones, which were analyzed 
above. Oblique dark solitons, which are located inside the outer Mach cone, they will be discussed in the next Section. 
One can also note that the oblique solitons are modulated by concentric circular waves, which were actually generated 
by the initial introduction (switching on) of the obstacle's potential. These circular waves are not related to the steady 
regime and will not be considered here. 
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We have compared analytical results ([24| for the crest lines with the numerical findings. The analytically predicted 
curves are shown by thin solid and dashed lines in Fig. 2. To simplify the pattern, we have chosen the parameters 
of the binary BEC so that the chemical potentials ([6]) of both components are equal to each other. As a result, the 
smaller sound velocity satisfies the relation = {gn — gi2)niQ — (522 — 912)^20 and hence the last term in Eq. (j25p 
as well as in analogous expression for n'2 vanishes. Therefore the ship waves in the region between the two Mach 
cones are not visible in agreement with the results of numerical simulations. Another linear mode describes the ship 
waves outside of the outer Mach cone, and the analytically calculated form of the crest line demonstrates excellent 
agreement with the numerically obtained one. 



IV. OBLIQUE DARK SOLITONS 



As mentioned above, in Fig. 2 one can observe oblique dark solitons located inside the outer Mach cone. It is worthy 
to note that they decay into vortices at end points, which is a result of the "snaking" instability of dark soliton stripes 
in 2D [13, m, H^l (see also [46*1). However, for large velocities U this instability is convective only [T^, which means 
that the dark solitons are effectively stable around the obstacle and, hence, the length of the dark solitons increases 
with time. The solitons originate from a depression in the density distribution formed behind the obstacle by the 
flow, therefore their depth also varies near the obstacle. However, upon sufficiently long evolution time, there exists a 
region where the oblique dark solitons can be considered as quasi-stationary structures. In this region, the solitons are 
described by the stationary solution of the GPEs. This allows us to take the stationary GPEs in the hydrodynamic 
form Q, 



V-(niUi) = 0, V-(n2U2) = 0, 
(Vtii)^ Ani~ 
8^^f 4^;""' (26) 



(ui • V)ui + 5iiVni + 512V712 + V ( ^ ^ - -j^ ) = 0, 



(U2 • V)u2 + gi2Vni + 522 Vn2 + V f - ^ ) = 0. 

\ 4n2 / 

Here, it is also taken into regard that the obstacle's potential is negligible far from it. Equations ([^5]) should be solved 
with the boundary conditions 

ni -> riio, n2 ^ riao, ui ^ (-t/,0), U2 ^ (-[/,0) at |x| ^ cx), (27) 

where —U is the common velocity of both components relative to the obstacle. Under the assumption that the solution 
depends only on 

X — ay 
Vl + a 



where a determines the slope of the oblique dark soliton, this system can be readily reduced to equations 

\{ni^(^ - 2nini_^^) + giin\ + gi2n\n2 + - {\q + y.i)n\ = 0, 



(29) 

\{'^2,(, - 2n2r7,2,^^) + 312^1^2 + 322"! + 5g"-20 - ihl + A*2)"2 = 0) 



where /ii and /i2 are the chemical potentials defined above in Eq. ([6]), and 



The flow velocities are related to the densities as follows: 

' ("-io + a^ni)U all (jiiQ — rii) 
(1 + a^)n^ ' (1 + a?)jM 



(30) 



i = l,2. (31) 



In general, system ([2^ has to be solved numerically. However, if the chemical potentials of the two components 
are equal. 



Ml = M2 = M, (32) 
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the system admits a simple analytical solution in a closed form. In this case, we look for the solution as ni = riio/(6)' 
^2 = ^20/(0 J reducing both equations to a single one, 

- 2//«) + f^f + - (ig + = 0. (33) 

Dark-soliton solutions to Eq. ([33|) are known [12]: 

1 - g/c^ 

ni = nis = nio/(0, "2 = = 7^2o/(0: /(6 = 1 r , -, , (34) 



cosh^ 



\Jc\- q {x - ay)/VTl~a? 



where we have also taken into regard that condition ([5^ leads to the following expressions for the sound velocities 
(for definiteness, we suppose here that 511,522 > 512, i-e., the inter-species repulsion is weaker than the repulsive 
self-interactions of the two components): 

= (511 - 512)^10, c+ = ^. (35) 

Obviously, solution ([34|) exists if the condition q < is satisfied. If we introduce angle 9 between the direction of 
the flow and the orientation of the dark-soliton stripe, so that a = cot 9, the latter condition can be transformed into 

sin^.<^^^, AU^E_. (36) 
Thus, the soliton must be located inside the outer Mach cone, which is defined by the equation 

sm9+ = ^. (37) 

Although we have arrived at this conclusion under assumption p2p . we conjecture that it is correct too in the general 
case of unequal chemical potentials, which is confirmed by the fact that our numerical simulations always produced 
oblique solitons confined inside the outer Mach cone. The respective numerically generated profiles of the densities 
are shown in Fig. [3] as a function of j/ at a fixed value of x; the corresponding phase profiles are also shown in the 
figure. It is observed that the oblique solitons are indeed located inside the outer Mach cone [as defined by Eq. (P7|) ]) 
but outside of the inner cone, which is defined by sin0_ = 1/M_. Profiles of the solitons' densities are close to the 
analytically predicted ones and the jumps of phases are also in agreement with the expected dark-soliton behavior. 

According to condition solution exists for any supersonic flow with t/ > c+, the same being true for the 
existence of numerical solutions to system (I29p in the general case, fj,i ^ ^2- However, that does not mean that the 
oblique dark solitons can be generated by any such flow. While the numerical results presented in Fig. [2] show that 
oblique solitons indeed exist for velocity U — 3.5, the opposite situation, when oblique dark solitons do not emerge, 
is presented in Fig. lUby means of density patterns which correspond to a lower supersonic flow velocity, U — 1.2. It 
is seen that vortex streets are generated in the latter case, rather than dark solitons. Such a behavior is related to 
the well-known instability of dark solitons with respect to transverse perturbations [13, [H, . The instability splits 
dark solitons into vortex-antivortex pairs, hence dark solitons cannot develop from the density depression behind the 
obstacle moving at a relatively low velocity. However, the numerical simulations presented in Fig. [5] indicate that 
the oblique solitons become effectively stable if the flow velocity is sufficiently high, as first was noticed in Ref. [l2| 
for the case of a one-component BEG. The stabilization was explained in Ref. [ij as a transition from the absolute 
instability of dark solitons to their convective instability, at some critical value of the fiow velocity, Ucr > c+, so that 
the unstable disturbances are carried away by the flow from the region around the obstacle where, as a result, the 
dark solitons look as effectively stable objects. Here, we aim to consider such a stabilization transition for the case of 
the two-component BEG, which features two unstable modes of perturbations around the dark soliton. 

To find the spectrum of small-amplitude linear waves propagating along the soliton, we now consider this solution 
in the reference frame in which the condensate has zero velocity far from the obstacle. To this end, we rotate the 
coordinate system by angle <p = arctana [recall a determines the orientation of the dark soliton, according to Eq. (p8|) ]. 
and perform the Galilean transformation to the frame moving relative to the obstacle at velocity {U cos(/?, U soup): 

X — xcos (p — y sin ip — U cos ip ■ t, , ^ 

■ ... (38) 

y — X sm ip + y cos ip — U sm ip ■ t. 

After the transformation, velocity fields (PT|) become 

ui = (w(nio/riis - 1),0) , U2 = (■L'(n2o/"-2s - 1), 0) , (39) 
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FIG. 3: (Color online.) Plots of densities (top) and phases (bottom) of the two BEC components for flow velocity U = 4.5. The 
positions of the Mach cones are shown by dashed lines, making it evident that the oblique dark solitons are located inside the 
outer Mach cone, in agreement with Eq. (I36p . The chemical potentials are ^ii = jj,2 ~ 3.5, cf. Eq. (|32p . and the configurations 
are shown at t = 30. 



and the densities take the form of 



nis = nio/(C), "-2s = n2of{C), 



/(C) = 1 - 



1 - v^/cl 



cosh 



— X — vt, 



(40) 



where the soliton's velocity in the new reference frame is 



U 



Below, we omit tildes attached to the new variables. 

We take small transverse perturbations of the dark-soliton solution as 

■01 = V'is(C) + (V'i + #1 ) exp {i(t>is{C) - ifJ-it) , 
= ^2s{0 + W + #2) exp {i(t>2s{C) - ifJ'2t) , 

where the unperturbed solution depends only on — x — vt, 



(41) 



(42) 



(43) 
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FIG. 4: Plots of densities of two BEG components for flow velocity U = 1.2 a.t t — 55. 



and phases cjjjs are related to the densities by equations 



Tljs 



(44) 



Perturbations V'' and i^" depend on y and t as exp{ipy + Tt). Substitution of expressions ((42|) into Eqs. ([T]) and the 
linearization with respect to ip' and itp" lead to an eigenvalue problem, 



/ Ai -Ln 

Lm Ai B 

^2 -L/2 

\ B Lr2 A2 



(45) 



A 



B = -2gi2^/nisn2s, 



" 2dC^ ~ Y~n^ '^2^'" ^P^^ ^^^"^^^ ~ ^'J"'" + '^J' (46) 

19^ 1 n,ow^ 1,2 2n 
= 2 ~ 2~^ ^ 2 - P ) - - gijnis + //j, 

J = 1,2, ;-i,2, 

System (|45|) determines the growth rates Ti,2{p) of small perturbations traveling along the dark-soliton's crest. The 
result is the presence of two unstable branches, examples of which are shown in Fig. [5l As one can see, both branches 
indeed feature regions of wave vector p with Rer(p) > 0. The transition to convective instability should be considered 
separately for each branch. 

Returning to the reference system attached to the moving obstacle, we get dispersion relations 

io^^2{p)^U,p + iT{p,v), (47) 



where Us — vsmip = aU j \J\ -I- is the component of the flow velocity along the dark soliton. The perturbations 
may be represented as Fourier integrals over linear modes obeying dispersion relations (j47p . 

/•OO 

<5ni,2 0c / <5ni,2(p)e*[™""'^'=^^^*lrfp. (48) 
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FIG. 5: (Color online.) Two branches of dispersion curves T{p) for unstable disturbances around the dark-soliton solution. Real 
parts of r are shown by crosses, and imaginary parts by dots. Parameters are gu — 1, (722 = 1-6, gi2 = 0.1, nio = 1, 7120 = 0.6. 



This representation implies that the integral is convergent, i.e., the wave packet is finite along coordinate y. Its time 
dependence at fixed value of y is determined by dispersion relations u;i,2(p)- Since expressions (|47|) contain imaginary 
parts, (5ni_2 can grow exponentially with time, which implies instability of the dark soliton, as is well known for 
the zero flow. Us = 0. However, it may happen that, for Us large enough, wave packets are carried away so fast 
that they cannot grow at fixed value of y, which is precisely the transition from absolute to convective instability 
[sot ■ Mathematically, the convective instability means that one can transform integrals over p in wave packets (j48p 
into integrals over w, because these wave packets have finite duration. In other words, function lo = w{p) can be 
inverted to define single- valued dependence p — p{lo)- Therefore, the distinction between absolute and convective 
instabilities depends on analytical properties of dispersion relations (|T7)) [sol . [3l|. Actually, the transition from 
absolute to convective instability is determined by critical points Pcr where dui / dp = 0, and function p — p{lo) changes 
its behavior: at Us < (C/s)cr is represented in the complex p plane by disconnected curves, whereas for Us > (C^s)cr 
these curves are connected with each other. In the latter case, one can deform the contour of the integration over p, 
with regard to the single- valuedness of w = (^{p), so as to transform it into an integral over uj. In other words, the 
spatial Fourier decomposition of the perturbation wave packet can be transformed to a temporal form, which means 
that the instability is convective. 

As is known, the asymptotic behavior of integrals (|48p is determined by branching points of function p = p{i-u), 
where duj/dp = 0. This yields equation 

Us = (49) 

where, as one can see in Fig. 5, T{p, v) has either real or purely imaginary values for real p. Therefore, critical values 
of Us, at which disconnected contours transform into connected ones, correspond to the appearance of a double root 
Pbr of Eq. on the real axis of p. This means dp^^/dUs = 00 &t Us ~ {Us)^^- The differentiation of Eq. with 
respect to Us then leads to equation 



d^T 
dp^ 



= (50) 

P=Pcr 



for the corresponding critical value, Pcr(w). The substitution of that value into Eq. (|^^ yields function Us{v). When 
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convective instability 
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FIG. 6: (Color online.) Curves separating regions of the absolute and convective instabilities of dark solitons in the two- 
component condensate with parameters gii = 1, 522 = 1-6, gi2 ~ 0.1, nio — 1, n2o = 0.6. The respective sound velocities are 
c+ = 1.03, c_ — 0.95. Squares and circles refer to two different unstable modes. Below each curve the corresponding mode 
is absolutely unstable and dark solitons cannot be created by the flow with velocity U less than Ucr — 1.5 for most slopes a. 
Dark solitons become just convectively unstable (effectively stable) above both curves, that is for U > Ucr = 1.5. 



this function is known, we find, with the help of relations 



the slope. 



and velocity. 



«=—=£==, Us = Usinip ^ ^E^= ^ av, (51) 
Vl+a vl+a 



ae,(«) = (52) 



U{v) = v^l + aliv), (53) 



as functions of v for all values in the interval (0 < u < c+). As a result, we obtain, in a parametric form, the 
dependence U{a) for the curve separating regions of the absolute and convective instabilities. Two such curves for 
both branches are shown in Fig. 6, where the region of the convective instability is located above both curves. It is 
seen that the dark solitons with large slopes a undergo the transition to the convective instability at flow velocities 
greater than Ucr ~ 1-5, as suggested also by the numerical findings reported above. 



V. CONCLUSION 



In this work, we have considered the effect of dragging a supercritical obstacle through a two-component Bose- 
Einstein condensate, which was motivated by a number of recent experiments in such settings, predominantly exam- 
ining the dynamics of mixtures of hyperfine states of ^^Rb. The presence of two speeds of sound in the mixture results 
in the existence of two Mach cones. The motion of the obstacle, in turn, produces two main features, namely, the 
linear "ship waves" and the oblique dark solitons. For the former pattern, we have developed a description of their 
density oscillations, which occur outside the Mach cones, in good agreement with numerical findings. On the other 
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hand, for the dark sohtons we have developed an analytical description of their profile, which was also found to be 
in good agreement with numerical observations. In particular, it was predicted that the dark solitons are confined to 
the area inside the outer Mach cone, which was confirmed by the simulations. 

This work may be a relevant starting point towards a more detailed understanding of the interplay between the 
inter-species interactions and the loss of superfluidity caused by the super-critical motion of defects. Among natural 
extensions of the present setting, there may be the consideration of the three-dimensional context, with vortex rings 
being formed as a result of the motion of the obstacle [which is closest to the experimental settings of reported in 
Refs. 0, [I]]- Another possibility would be to consider the supercritical motion of the obstacle in a three-component 
spinor mixture (47| , where it would be relevant to examine the role of the spin-dependent and spin- independent parts 
of the interatomic interactions in producing patterns such as those considered herein. Such studies are currently in 
progress and will be reported in future works. 
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